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Abstract. Various classification schemes exist for homo- 
geneous and isotropic (CP) world models, which include 
pressureless matter (so-called dust) and Einstein's cos- 
mological constant A. We here classify the solutions of 
more general world models consisting of up to four non- 
interacting fluids, each with pressure Vj, energy density ej 
and an equation of state Vj — (jj — \)tj with < jj < 2. 

In addition to repulsive fluids with negative pressure 
and positive energy density, which generalize the classical 
repulsive (positive) A component, we consider fluids with 
negative energy density as well. The latter generalize a 
negative A component. This renders possible new types of 
models that do not occur among the classical classifica- 
tions of world models. Singularity-free periodic solutions 
as well as further 'hill-type', 'hollow-type' and 'shifting- 
type' models are feasible. 

However, if one only allows for three components (dust, 
radiation and one repulsive component) in a spatially flat 
universe the repulsive classical A fluid (with A > 0) tends 
to yield the smoothest fits of the Supernova la data from 
Perlmutter et al. (1999). Adopting the SN la constraints, 
exotic negative energy density components can be fittingly 
included only if the universe consists of four or more fluids. 
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1. Introduction 



While in the Newtonian description gravity is only con- 
sidered as an attractive force, repulsive effects are pos- 
sible in general relativity, as e.g. caused by Einstein's 
(1917) cosmological constant A (if positive). The variety 
of general-relativistic CP world models (CP=cosmological 
principle which requires homogeneity and isotropy after 
averaging on suitable scales; we require exact homogeneity 
and isotropy for the mathematical CP models) including 
the A-term and non-relativistic incoherent, comoving (i.e. 
with pressure V — 0) matter (so-called dust) was classi- 
fied using different parametrizations by Robertson (1933), 
Stabell & Refsdal (1966) and Carroll et al. (1992). From 



the seventies till the nineties, a widely followed assump- 
tion was to set A = 0, partly promoted by the inability of 
current quantum field theory to predict a realistic value 
of A (e.g. Weinberg 1989). 

This attitude changed during the last few years. The 
near-flatness suggested by cosmic microwave background 
(CMBR) observations (Lineweaver 1998, de Bernardis et 
al. 2000) and the time-scale problem caused by a tendency 
to a large Hubble constant led an increasing number of 
workers to consider world models with A or similarly re- 
pulsive fluids (e.g. Ostriker & Steinhardt 1996; Turner & 
White 1997). Recently, direct motivation for such an ap- 
proach has come from observations of supernovae of type 
la, taken as cosmic 'standard candles', that require the 
presence of an accelerating component in the expanding 
universe (Riess et al. 1998; Perlmutter et al. 1999), primar- 
ily taken as evidence for the action of A. If so, CMBR data 
and peculiar motions of galaxies can in principle break 
the degeneracy between the decelerating and accelerating 
components given by the SN data alone, so that it may be- 
come feasible to determine an accurate value of A (Zehavi 
& Dekel 1999). 

However, given our ignorance about the physical ori- 
gin of A it may be premature to identify the A term with 
the repulsive component. E.g., in the early universe in- 
flaton fields or domain walls are well known to serve as 
expansion-driving forces as well (Kolb & Turner 1990). 
As will be recalled in Sect.||, any ideal fluid with negative 
pressure may force expansion of space if its equation of 
state obeys a certain condition. 

It is the goal of the present work to extend the above 
cited classifications of CP world models in two ways. 
Firstly, we allow for more than two components. Secondly, 
we allow for a range of fluids which have in common that 
each equation of state yields proportionality between en- 
ergy density e and pressure V. This form of the equation 
of state includes the classical cosmological matter compo- 
nents dust, radiation and the A fluid, but also encompasses 
more general repulsive and attractive components. In ad- 
dition to the above mentioned negative-pressure compo- 
nents (generalizing the A > case) we also allow for fluid 
components with negative energy densities which general- 
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ize the case of a negative A. The present work extends a 
discussion started in Schulz (1998). 

2. Basic equations 

To fix our notation we now compile the basic equations 
for the subset of CP world models to be studied here. The 
vacuum velocity of light is set c = 1. Firstly, we recall that 
the A term is equivalent to a fluid with constant energy 
density 6a. 

2.1. Constant- density fluid 



Einstein's field equation in the general form reads 



(i) 



in which we take as source term on the right side a sum 
over ideal-fluid energy-momentum tensors of the form 



Vg- 



(2) 



V and e are pressure and energy density measured in a lo- 
cal comoving frame, in which the four-velocity reduces to 
vP = 1, v?~ = u 2 = = 0. Despite lacking a definition of 
their four-velocity, isotropic relativistic components (null- 
fluids) can be subsumed under this scheme with an equa- 
tion of state V = e/3. The A term in Eq.(|j]) can be omitted 
on the left side and instead, on the right side, be repre- 
sented by a 'fluid' after substituting A = 87rGeA with a 
constant energy density e\ and adopting an equation of 
state Va — — £a for the A-fluid. 

2.2. Fluid equations 

The basic symmetries given by the assumption of homo- 
geneity and isotropy yield the Robertson- Walker (RW) 
line element which reads with spatial polar coordinates 

-'■ - <7, x 2 = 8, X 3 = <f) 



dr 2 = g^dx^dx" 
= dt 2 - a 2 {t) 



da 2 



1-ka 2 



+ a 2 (d9 2 + sin 2 6>d</> 2 ) (3) 



k € {0, — 1, 1} is a curvature index. As is common practice 
we call models with k = 0, —1, 1 flat, open or closed, re- 
spectively, without emphasizing its relation to spatial cur- 
vature. Here a is fixed for comoving fluid elements. The 
right side of Eq. (|l|) is subject to the conservation law for 
matter fields T 11 ".^ = yielding the equation of continuity 
for perfect fluids 

which reduces with the RW line element to the cosmolog- 
ical energy equation 



e=-3(e + 7>)^ or ^l = -3Va 2 . 
a da 



(4) 



The relativistic Euler equation (also stemming from 
= 0) 



{e + P)u' l u v . Vk = [ a r-vf l u v )V., ll 



(5) 



reduces to u tJj u' J -^ = 0, showing that an ideally homoge- 
neous and isotropic universe has geodesic flow lines. 

The relativistic analogue of Newtonian mass density 
in the Euler equation Eq. (||) is (e + V) and can therefore 
be considered as generalized inertia. For the stability of 
'normal' matter (e + V) > is required. Fluctuations in 
the spatial distribution of energy cause deviations from 
geodesic flow, if (e + V) ^ 0. For the A-fluid, because 
of its vanishing inertia, inhomogencitics arc inconsistent 
with Eq. (|). 

Inserting the RW line element Eq. (||) into the field 
equation Eq. (hi) leads to the Friedmann equations 



a 
a 

and 



4ttG 



(e + 3V) 



(6) 



(7) 



H(t) is the Hubble parameter at time t. The system of 
Eqs. (Q), (Q) and (Q) is redundant; from any two of these 
equations one can derive the third. 

2.3. Equation of state 

For full specification of a gravitating source one needs an 
equation of state. Here we assume for each cosmic fluid 



(7 - l)e 



(8) 



where 7 is a constant specifying the fluid. This form was 
chosen by Zeldovitch (1962) as an equation of state also 
common for nuclear matter. It applies to the best known 
cosmic fluids like (e.g. Kolb & Turner 1990) 

— dust (V = 0; 7 = 1), 

— radiation and any other extremely relativistic matter 
{V = e/3; 7 =4/3), 

— the constant-density A-fluid (V = — e;j = 0), 

— a network of slow cosmic strings (V = —e/3; 7 = 2/3) 
or 

— domain walls (V = -(2/3)e; 7 = 1/3). 

Apart from incoherent matter (dust) and radiation the 
physical background of cosmic fluids with equation of 
state (||) is hardly understood. Therefore, in case of 7 7^ 1 
or 4/3 we use the designation exotic fluid. 

We confine the parameter of the equation of state to 
the range 



< 7 < 2 implying c s < c 



(9) 



where the speed of sound c s is restricted not to surpass 
the velocity of light c. 
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Inserting Eq. (g) into Eq. (0) and integrating yields 

-3 7 



e = e — 



(10) 



where the constants eg and a are here fixed at t = i , the 
present epoch. 

For the A-fluid (e A + V A ) = Eq. (g) yields e A = or 
€\ = const as required in Sect. 2.1. Eq. (^|) rules out any 
pressure gradient in the A-fluid. Hence, the density of a 
universe consisting solely of the A-fluid would be spatially 
and temporarily constant, i.e. it would satisfy the Perfect 
Cosmological Principle of Bondi & Gold (1948) and Hoyle 
(1948). Indeed, McCrea (1951) derived the equation of 
state V = — e for the matter content of the steady-state 
universe. 

Since the definition of the cosmological constant by 
Einstein (1917) positive as well as negative values of A 
have been considered. Transferring this to the fluid de- 
scription given in Sect. 2.1 suggests to allow for positive 
and negative energy densities of other 7 fluids in our phe- 
nomenological approach as well. 

Static world models like Einstein's (1917) historic so- 
lution will not be included in th e cla ssification diagrams 
below (albeit mentioned in Sect. ^2). The FLRW inter- 
pretation of the redshifts of galaxies makes such models 
obsolete. In our context, de Sitters (1917) historic world 
model corresponds to a single-fluid universe with 7 = 0. 
Despite misleading information in some textbooks, Fried- 
mann's (1922, 1924) solutions are two-component models 
consisting of dust (7 = 1) and the A fluid (7 = 0) (as has 
been noted by Felten & Isaacman (1986) as well). 

2.4- The multifluid approach 

We assume a finite number of non-interacting] 1 ! cosmic flu- 
ids with equations of state of the form Eq. (|T, an energy 
equation of the form (g) and, hence, Eq. (10|) holds for 



each fluid so that 
e = ^2 ei V = y^(7,— l)e< with e. L = e oi 



-3 7 i 



(11) 



Defining the normalized density f2i, the normalized scale 
factor x(t) and the parameter of state cti by the equations 

SttG . . aft) , „ „ 

:= -r-^rreoi, x(t) := and a; := 2 - S'fi 



3HS 



«o 



we can write Eqs. (g) , (0) in the form 
ctiVliX a% =: fi(x) 



My, 

2x ^ 



x 







n, 



X 



(12) 



(13) 



1 Non-gravitational interactions, i.e. interactions other then 
through the Friedmann equations, are excluded. 



In these equations we have introduced the abbreviations 
/1 (x) and ji (x) for the right-hand sides of the Friedmann- 
equations. Therebye we emphasize that f\ and fi are ex- 
plicit functions of the scale factor x, which only mediate 
the dependence on cosmic time t. It is mathematically 
crucial for the classification that f\(x) and f2{x) do not 
explicitely depend on time. We also want to distinguish 
between the different role of the dynamical functions x 
and x and the conditions fi(x) and f 2 (x) that determine 
the temporary behavior of the scale factor x(t). 

To derive Eq. ( |l3|) one has to use the present-epoch 
(t = t ) version of Eq. (0) to eliminate the curvature index 
k of Eq. §) 



a Q H Q 



5> 



(14) 



By differentiating Eq. ( |13[ ) with respect to x one ob- 
tains 

ld/ 2 (x) ldx 2 (x) 



2 dx 2 dx 



x(x) = fi{x), 



which is an important property to be used below. 
Eq. (jll|) can be formally integrated by 



Ho(t — tinit) 



dx 



(15) 



(16) 



In big-bang models the initial time ti n a = can be chosen 
for the 'big bang' (given by the first instant with x = 
when progressing backwards in time). In models with an 
infinite past t; n it can be chosen in accordance with one's 
wish how far to follow the model. Eq. ([!(]) gives t(x) which 
by inversion (sometimes of several monotonic branches) 
leads to the desired world model x(t). 

In practice, rather than carrying out the direct integra- 
tion it is frequently more convenient to solve Eq. (|^) via a 
standard fourth-order Runge-Kutta routine. Good exam- 
ples for the computation of the case A + dust models are 
given by Felten & Isaacman (1986). 

3. Cosmic dynamics 

In the following we call components with a > A-like, 
those with a < dustlike because of their dynamical sim- 
ilarity to the A-fluid (a = 2) and cosmic dust (a = — 1), 
respectively. 

3.1. Boundaries of the scale factor 
The physical meaning of x(t) requires 



if > 



(17) 



restricting x 1 to the positive branch of f2(x), so that x 
and x are real valued. We assume x > and identify 
the big-bang with x = 0. As a consequence of Eqs. @ 
and (||), fi and f 2 (cf. Eqs. ([l2]), (|l^)) depend on x and 



4 



J. Thomas & H. Schulz: Multifluid world models 



only indirectly on cosmic time t via x(t) (as pointed out 
in Sect. 2.4). Thus, geometric features of the solution x{t) 
are naturally related to x rather than to t. For instance, 
if at a certain time i max the cosmic evolution passes a 
maximum x max = x(t max ), then for any other maximum 
at t' > t max necessarily x(t') — x maK holds, meaning that 
£max is a global and unique maximal value of x(t). Obvi- 
ously the same holds for a minimum x min . Consequently, 
there is at most one maximal value of x or at best only 
one minimum, which we denote by x max and x min , respec- 
tively. It is nevertheless possible that x(t) obtains these 
values infinite times (in periodic models). 

For the classification of solutions it is appropriate to 
distinguish between two different types of maxima or min- 
ima. Firstly, let us consider a maximum of the scale factor 
which is a simple (or first-order) zero a; max > 1 of fa(x) 
by requiring (cf. Eq. (|l5|)) 



f2{x m nx) = 0, fi(x maK ) < 0. 



We choose an initial instant of time ij with X{ = x(U), 
at which the universe is assumed to expand. During sub- 
sequent evolution the universe reaches the maximum at 
time 



dx 



(19) 



Extracting the simple zero of f 2 {x) we define a function 
g{x) by 



f 2 (x) = (x max - x)g(x), 

so that according to Eqs. (p7|), ( |l8| ) 

g(x) > for Xi < x < x n 



(20) 



Thus, 



is finite and 



^m aX ti 



c := sup 



1 



dx f dx 

< c 



\fx max x\J 'g(x) 



V^max X 



— 2c-\/x ma x Xi- 

Hence, in this case x max will be reached after a finite 
duration of time and may therefore be designated as a 
maximum (in contrast to the asymptotic upper boundary 
described below). After this finite i max the universe has 
to contract. World models may exist which pass through 
such a maximum infinite times. 
Analogously x m i n < 1 with 

/2(Zmin) = 0, /l (a: m in) > 

will lead to a minimum in the above sense (that could 
occur before i,-). 



A further possibility is a value x up > 1 with the prop- 
erties 



f2(x up ) = and fx(x up ) = 0. 



(21) 



Again, from an initial scale factor Xi = x(ti), where the 
universe is assumed to expand, it will reach a; up at t up 
in accordance with Eq. (|l9|). However, when following the 
above steps and defining a function g{x) such that 



inf 



1 



x£[xi,x up ] \J g{x) 
is finite, we now obtain 

f 2 (x) = (x up -x) n g(x), (22) 
with n > 2. Together with Eq. ( ^9|) this yields 



(18) ^up w 



dx 



> c 



(x up - x) 2 y/g(x) 

P 

dx 

(^up 



> c 



dx 



(23) 



j up 



where the last inequality holds at least in an appropriate 
neighborhood of x up . The last integral in ( p3|) diverges 
and, consequently, x np will be reached only asymptotically 
after an infinite time, meaning t up — » oo. 

Furthermore, the scale factor approaches a quasistatic 
limit at a; up . To show this we apply Leibniz' rule on 

dx dx . 
dt dx 



With := one obtains 



d n+1 x (dx ± 



(n) 



dt n+1 \dx' 
which can be written as 



k=0 



(n-fc) 



(X) 



(fc) 



Jn+3) 



^ \k 

k=0 



(X) 



(k) 



(24) 



So, with vanishing 

f2(x up ) = and fi(x up ) = 

all higher time derivatives of x(t) vanish at x up . (Note that 
x and its derivatives are continuous for x > 0.) 

In the case corresponding to Eq. ( pi]) we call a; up an 
asymptotic upper boundary. Analogously if xi ow < 1 and 

/2(Tlow) = 0, /i(xi ow ) = 

then we call cci ow an asymptotic lower boundary. The uni- 
verse can begin or end with such states. 

It is worth mentioning that due to Eq. ([l5]) a maximum 
(or a minimum) of x{t) is obtained where fz{x) has a 
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simple zero. An asymptotic boundary is obtained by a 
higher-order zero. 

Since the distinction is clear from the context we use 
below the designations x m i n for x m i n and x\ ow and x max 
for x max and x U p. 

The interval within which x(t) is bounded for all pos- 
sible times t is the interval [x m in , i£ max ] that comprises 
unity and is limited by the first zero < x m { n < 1 of 
fi(x) {x m i n :— 0, if no such zero exists) and the first zero 
£max > 1, respectively (x max := oo, if no such zero exists). 

3.2. Basic solution types 

We first assume a currently expanding universe, i.e. Hq > 
(as the real universe appears to be). Depending on 
whether has zeros of the above described types or not, 
we can distinguish nine basic types of worldmodels: 

1. There is no zero of fs(x) with < x m - ul < 1. 

Such a universe starts with a big bang because, when 
progressing backwards in time, its slope x stays posi- 
tive so that eventually x = will be reached. If there 
exists 

(a) no zero x max > 1 of fs(x), then the universe ex- 
pands forever and the scale factor grows to infinity; 

(b) a higher-order zero X max > 1 of /^(s), then the 
expansion will stop after an infinite time and the 
universe becomes asymptotically static with x — ► 

(c) a simple zero a; max > 1 of f2{%), then x(t) will rec- 
ollapse to x = after passing the maximum cc max . 

2. fi{x) has a higher-order zero < x m i n < 1. 

In this case the universe starts with an asymptotically 



static state at x 



If there exists 



(a) no zero x max > 1 of f%(x), then the universe ex- 
pands forever and the scale factor grows to infinity; 

(b) a higher-order zero X max > 1 of f2(x), then the 
universe becomes asymptotically static after infi- 
nite time with the asymtotic scale value x — ► x max ; 

(c) a simple zero x max > 1 of f2(x), then x(t) reaches 
a maximum and after that decreases. The contrac- 
tion stops after an infinite time with x — > Xmin and 
the universe becomes asymptotically static. 

All these models would have begun in the infinite past. 
3. There is a simple zero < x m i n < 1. 

The universe has not started with a big bang. If there 
exists 

(a) no zero x max > 1 of f2(x), then the universe con- 
tracts from infinity, passes a minimum and expands 
back to infinity: 

(b) a higher-order zero a; max > 1 of f^ix), then the 
universe starts with the asymptotically static state 
x = £ m ax, contracts to a minimal value x m i n , ex- 
pands again and after infinite times reaches the 
asymptotic value £ max ; 



(c) a simple zero x max > 1 of f2(x), then x(t) is peri- 
odic. 

Because there is no big bang beginning, such world 
models are infinitely old. 

Solutions with at most one boundary are well known 
from the classical dust models with Einsteins cosmological 
constant A and can be found in Stabell & Refsdal (1966). 

Models trapped between two boundaries are shown in 
Fig. [I]. There are 'shifting' (pb| ), 'hill- type' (pc|), 'hollow- 
type' ( |3b| ) and 'oscillating' ( pc] ) solutions. Apart from the 
periodic type they require more than two fluid components 
and were not described in earlier classifications. 

The case Hq = leads to one more solution class: static 
universes with x — for all times. One famous member 
of this class is Einstein's static universe with dust and 
cosmological constant. 

If Hq < 0, then the classification remains qualitatively 
the same as for Hq > 0. However, solutions with Hq > 
0, that have a symmetry-point t sym with x(t sym + t) = 
x(t S ym — t) will be translated in time, when Hq becomes 
negative. This is the case for the original classes 'big bang 
+ max' ([lc|), 'hill- type ' J2c| ), 'bouncing' (|3^), 'hollow-type' 
( |3b| ) and 'oscillating' (|3c|). 

On the other hand, those solutions for which x(x) > 
holds for all x e [x m i n , x max ] are reflexions of solutions 
with Hq > with respect to the (t = £o)-axis. This 
concerns the original classes 'big bang without recol- 
lapse' (|la|), 'big bang + asymptotic upper boundary' (|lb|) , 
'asymptotic lower boundary without recollapse' ( pa|) and 
'shifting' ©. 

3.3. Repulsive and attractive fluids 

Eq. ( |T2] ) shows that a fluid is attractive, as long as aVl < 0, 
repulsive if aVl > 0. If fi > 0, A-like fluids are repulsive, 
dustlike are attractive, whereas with fl < A-like fluids 
are attractive and dustlike are repulsive. 

However, while attractive dustlike components pro- 
mote a x — singularity, attractive A-like fluids tend to 
limit the expansion to an upper boundary x max . This can 
be seen from Eq. ([l3]) in which each fluid, in addition to 
affecting the curvature with a HqVI term, adds a term of 
the form HqHx 01 to the cosmic scale factor. For dustlike 
attractive fluids a < and the term is positive. Moreover, 
the component enforces contraction when the scale factor 
is decreasing. Conversely, for A-like components a > 
but for attraction < 0. The term is negative and its 
influence increases, when the universe expands. 

Repulsive dustlikc fluids counteract the formation of a 
singularity (it depends on the effective sum over all fluids 
whether a singularity will be truly avoided), repulsive A- 
like work against formation of an upper boundary of the 
scale factor in an expanding universe. 

For a dustlike component Qx a vanishes when x — > oo, 
and the component only contributes to curvature. The 
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x(t) 1- 



x(t) 1 




Fig. 1. When fluids with negative energy are present x(t) 
can take one of the above forms. Oscillating solutions 
are possible with two components, while the shifting-type 
(upper part, two asymptotic boundaries) as well as the 
hollow- and hill- type (lower part, one asymptotic bound- 
ary) require more than two fluids. 

same holds for A-like fluids when x — ► 0. Therefore world- 
models, that consist of fluids that are all dustlike or all 
A-like are dominated by curvature effects. 

Models with exclusively dustlike, attractive fluids have 
necessarily a big bang. If the geometry of the model is 

— open, then the universe expands forever ( |la| ) 

— flat, then the universe expands forever and becomes 
asymptotically static with x — > oo (lb) 

— closed, then the universe passes a maximum and col- 
lapses to x — Jlc|). 

Worldmodels that contain only dustlike, repulsive fluids 
are always open and of the bouncing type ([3a]). 

Analogously models with exclusively A-like, repulsive 
fluids expand forever, where 

— open models contain a big bang (|la|), 

— flat models contain an asymptotic big bang (|2a|) and 



— closed models have a minimum of the class (|3a|) . 

Universes consisting of only A-like, attractive fluids are 
necessarily open and from the type (|lcj). 



4. Classification of two-component models 

For classifying two-component models in the allowed range 
of fluid parameters, one needs to find the zeros of f2{x). 
For models with two cosmic fluids (fii, ct\) and (r2 2 ) <2 2 ) 
the solution types occupy well-defined areas in the (f2i, 
f2 2 )-plane for prechosen ct\ and a 2 - Roots x c of the equa- 
tions fi{x c ) — f2(x c ) — lead to the boundary curves 
that separate different types of solutions. 

In the following, the parameters of state ai and a 2 



as defined in Sect. 2.4 are regarded as arbitrary constants. 
Without loss of generality we assume > |a:2 1 - 

4-.1. Transition curves 
Eqs.@, @ reduce to 



Hi (Q lX ai + Q 2 x a2 + 1 - O x - Q 2 ) 



tt2 

x = -^- (ai^ix Ql + a 2 fl 2 x a2 ) . 
2x 



(25) 
(26) 



Solving these equations for fi(x c ) = /^(^c) = yields 
n 2 



a 2 (l 



L )-«i( 



n 2 = 



Oil{X c 

01 



ai(xc' 2 - 1) + 0:2(2;° 



(27) 



(28) 



Eqs. (|2~7j), (28) define a non-continuous curve in the 
(fli , O2) plane parameterized by the critical scale factors 
x c . The discontinuity at x c = 1 divides the curve into two 
parts tracing models with an asymptotic upper boundary 
[x c > 1) or an asymptotic lower boundary (x c < 1). Each 
branch of the curve marks a transition between regions 
with different solution types. To show this we study the 
case of two dustlike fluids in detail. Other cases can be 
treated in an analogous way. 

4-2. Two dustlike components 

Recall that a dustlike fluid is attractive, as long as f2 > 0. 
The solution given by a chosen pair (fli, f^) depends 
strongly on which component is attractive, which is re- 
pulsive. 

ill > 0, > 0: / 2 has no extremum and the type of 
solution depends only on the curvature. We find 



lim f 2 = Hi (1 - fix - fi 2 ) 



(29) 



and conclude that closed models (k > 0, cf. Eq. (|14|) ) nec- 
essarily lead to a zero of f%{x) greater than unity. Thus 
they have a maximum and recollapse into a big crunch. 
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Open and flat ones expand forever, the latter stopping 
their expansion asymptotically. All models start with a 
big bang singularity. 

Hi < 0,^2 < 0: All models are open (fc < 0, cf. 
Eq. ([lj)) and have no maximum because of Eq. (^9|) . Af- 
ter collapsing from infinity and passing over a minimum 
they expand back to infinity, because 

lim f 2 = -oo if Qi < (30) 

shows that there is a zero of f 2 (x) in the interval [0, 1]. 

£li < 0, f^2 > 0: From Eq. ( |30"| ) we see that such models 
cannot collapse to x = 0. /i has only a single zero and 
therefore f 2 one extremum. According to Eqs. (|2^), ( |30| ) 
it must be a maximum, what geometrically means that 
f 2 has no, one or two zeros, depending on whether f 2 is 
negative, zero or positive at its maximum. There is at 
least one zero, because f 2 (to) > excludes a negative 
maximum. Summarizing, if the universe is open or flat it 
is of the bouncing type (|aj) (Eq. (|9|) , only one zero) , but 
in closed worlds the solution is periodic. 

ill > 0, < 0: Again f\(x) has one zero, but now 
leading to a minimum of f 2 (x), because of 

lim f 2 = oo if O x > (31) 

x^Q 

and Eq. (||^). Therefore f 2 {x) has one zero greater unity if 
Qi+Q 2 < 1 and the solutions begin with a big bang, pass 
through a maximum and recollapse. If Qi + CI2 > 1 then 
f 2 {x) has no, one or two zeros, and now all possibilities are 
allowed in the sense of Eq. @. If there is no zero, then 
the universe starts with a big bang and expands forever. 
Two zeros lead to solutions with one extremum, because 
either both are less or both are greater than unity. Finally, 
if f 2 {x) has only one zero, then this must be the zero of /i 
and, hence, fii, lie on the curve defined by Eqs. (p7|), 
(p8|). Thus, crossing this curve the solutions will skip from 
one type to the other. 

Fig. U summarizes the results for specific values ct\ — 
— l.7a 2 = —1. If ai becomes smaller, or a 2 becomes 
larger, then the transition curves become steeper. In this 
sense, one may refer to a as a parameter characterizing the 
strength of the attraction or repulsion of a fluid. The qual- 
itative properties of the transition curves do not change 
under variation of ot\ and a 2 . 

4-3. The general case 

The general case of two arbitrary components follows anal- 
ogously. 

Fig. U shows the regions of solution types for two A-like 
fluids. The assumption ct\ > a 2 enforces the transition 
curve to occur in the lower right quadrant. The curves 
look similar to those of Fig. pi however the roles of positive 
and negative energy densities have interchanged so that 
the domains 'bouncing' and 'big bang+max' flipped over. 



^ ^ oscillating 

^ 2- 
bouncing ^ ^ 

-2 


big bang + max 
- . 2 fl l 




-2- 

bouncing 

-4- 
-6 


\ big bang + max 

\ big bang 

\ without recollapse 

bouncing \ 



Fig. 2. The density-parameter plane for two dustlike fluids 
with (fix, cti = —1.7), (f^2, a 2 = —1). The type of solution 
depends on the strength of each component. Models on the 
dashed line are flat. 



^ oscillating 
big bang + max v ^ 

* 


a 2 

bouncing 

2 


-2- 

big bang + max 

-4 
-6 


\ bouncing 

\ big bang 

\ without recollapse 

- big bang + max \ 



Fig. 3. The density-parameter plane for two A-like fluids, 
(fii, ai — 1.7), (fl 2 , a 2 = 1). The roles of positive and 
negative energy densities compared to Fig. || have inter- 
changed. As in Fig. [2] flat models are related to the dashed 
line. 



This illustrates the opposite behavior of dustlike and A- 
like fluids. The transition curves will get steeper, when 
either ot\ becomes larger, or a 2 becomes smaller (or both 
together). 

In Fig.|] a dustlike fluid (x-axis) is mixed with a A-like. 
Now the transition curves have moved into the upper right 
quadrant. This reflects the fact that fluids with different 
signs in the parameter of state act contrary when the signs 
of the densities are equal in contrast to the previous cases. 
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/ big bang without recollapse 

^^^^^^^^ 


-2 
oscillating 

-2 
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Fig. 4. The case of one dustlike fluid (Sli, a\ — —1.4) 
and one A-like fluid (f^2 , «2 = 1-5) shown here is similar 
to that of universes containing dust and a A-term. The 
transition curves become steeper if a\ or a.% decreases. 



a 

big bang + max 

1 

-l f 5 


\ bouncing 

big bang without recollapse _______ 

, M ? " 


big bang without 
■ — — ___^ recollapse 

^\^1 
-2- 

bouncing 

-3 


big bang 

without big bang + max 
ecollapstf 



Fig. 5. Possible solution types for three components, 
among which dust (_-m = 0.3) and radiation (__ 7 = 
5.9 10 -5 ) are fixed. SI and a are the parameters of a third 
component. Models on the dashed line are flat. 



5. Models with one exotic component 

Recent observations of supernovae of type la (Perlmut- 
ter et al. 1998) suggest the presence of at least one re- 
pulsive fluid in addition to the cosmic radiation back- 
ground (CMBR) (fly , a 7 = —2) and pressureless (inco- 
herent) matter (dust) (SImj^m = — !)• For f_ 7 we can as- 
sume Sl 7 £ [0,1]. Incoherent matter encompasses 'ordi- 
nary' baryonic matter and possibly what is called cold 
dark matter (see Sect.|| for details). Due to our ignorance 
of the physics and amount of dark matter it seems appro- 
priate for classification not to restrict the value of __m- The 
incoherent component is attractive, as long as its density 
is positive. 

For three components Eqs. (|l2|), ( |l3| ) read 







Hi ( — + "-=!■ + n x a + 1 

x x l 



SL 



n 



M 



si 7 - n 



H 2 



x 



20, 



= [ --^-^L+aSlx" 
2x 



(32) 



(33) 



As in Sect.jH we need to find out simultaneous zeros of 
/2 and /i for classifying the solutions. However, we now 
regard the parameters of the 'third' component, param- 
eter of state a and normalized density Q as independent 
variables. The corresponding equations for (SI, a) cannot 
be solved analytically for a. 

By transforming Eqs. (|32|), ( |33"| ) into 

^m/x c + 2S1 7 /^ 



Q(a,x c ) = 



(34) 



ax 



a+2 



(QmX c 



X, 



20 7 ) 

n 



M 



si 7 = o 



(35) 



one can proceed to solve Eq. (|35| ) numerically for a(x c ). 
Inserting this result into Eq. ( |34[ ) yields Sl(q(x c ), x c ). 

If __m > and S1m + Sl 7 < 1, then Eq. ( p5[ ) has no 
solution for x c > 1 and the curve £l(x c ) = £L{a(x c ),x c ) 
describes models with an asymptotic lower boundary, sep- 
arating bouncing solutions from those with big bang, that 
have no maximum. It splits into two branches for a > 
and a < 0, respectively. Fig.|| shows the distribution of 
solution classes in the (SI, a) plane for this case. The indi- 
vidual domains arise from arguments comparable to that 
in Sect. 4.2. Due to the restrictions we made for S1m and 
fi 7 none of the new solution types of Fig.[l] arises here. 

In Fig. ^ the result of classification is shown for S1m > 
1 — S7 7 . The greater value of S1m increases the number of 
posssible solution types. By the differences of Figs.|| and 
H the general fact is expressed, that a set of dustlike fluids 
fli > obeying ]T\ fij < 1 cannot lead to a maximum 
of x(t) without further attractive components. Therefore 
Fig. |6| shows solutions with a maximum, where the third 
component is repulsive (upper right quadrant), that are 
caused by the greater value of __m- Accordingly, there are 
two different types of transition curves (as in Figs . §-Q) • 
Between 'bouncing' and 'big bang without recollapse' the 
solutions have an asymptotic lower boundary, between 'big 
bang without recollapse' and 'big bang + max' they have 
an asymptotic upper boundary. 

Fig. [7] shows the solution types, when S1m < Sl cr i t , 
where Sl cr j t represents the critical dust density that ex- 
cludes a big bang in dust-radiation-models. Sl C rit depends 
on Sl 7 , for example with Sl 7 = 5.9 10~ 5 (see Sect.^) one 
finds fi c ri t = -0.015. 

For the transition between different solution types the 
same holds as for Fig.pl. 



J. Thomas & H. Schulz: Multifluid world models 







big bang - 
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bouncing oscillating 
-3 



big bang 

without recollapse 




big bang + max 



bouncing 



bouncing 



Fig. 6. Classification of solution types for three compo- 
nents, f2 7 = 5.9 10~ 5 , f^M > 1 — and an arbitrary third 
component (f2, a). The models on the dashed line are flat, 
that on the bold part of the transition curve have hill-type 
solutions. 



oscillating 



-1- 



bouncing 



-3- 



bouncing 



bouncing 



big bang 
without 
recollapse 



oscillating 



big bang + max 



1-5 



Fig. 7. Solutions with three components, £1 7 = 5.9 10" 

< —0.015 and no restrictions for the third component 
(f2, a), flat models are on the dashed line, hollow-type 
solutions on the bold part of the transition curve. 



6. Models with two exotic components 

An increasing number of components makes the classifica- 
tion more and more difficult. When the geometry of f2(%) 
becomes richer in variants, e.g. if there are more extrema, 
then the method of classification introduced above is no 
longer useful, the simultaneous zeros of / 2 (x) and fi(x) 
are in general no longer distinctive for different solution 
types. However, one can get the solution type directly from 
/2(x) in any case (see Sect.jsJ). 

As in Sect.|| we represent the solution types in the 
(fii, 02)-plane of the densities of the additional compo- 



Fig. 8. The possible solutions in four component models 
with dust, radiation and two dustlike, exotic components 
(fii, ai = —0.9), (0 2 , ax = —0.5). Solutions on the tran- 
sition curves have an asymptotic lower boundary, on the 
bold part they have also a maximum (hill- type). Models 
on the dashed line are flat. 



ncnts ax), (O2, 0:2)- They now look like in Figs.|q-|lO|, 
for which we have adopted f^M = 0.3 and f2 7 = 5.9 10 -5 
(see Sect. [7]). However, there is no qualitative change in 
the classification as long as Dm + 7 < 1 and both densi- 
ties are positive. Each of the figures corresponds to one of 
Figs. ||-[| according to the combination of dustlike and A- 
like fluids as additional components. As a difference, the 
transition curves are no longer restricted to one quadrant 
(though this can happen with suitable values of ax and 
CK2). The bold parts of the curves mark the border be- 
tween oscillating solutions and those with big bang and 
maximum. They represent 'hill-type' solutions. The tran- 
sition curves between 'bouncing' and 'big bang without 
recollapse' belong to solutions with an asymptotic lower 
boundary, those between 'big bang without recollapse' and 
'big bang + max' belong to solutions with an asymptotic 
upper boundary. 

With increasing J7m or 7 the regions of big-bang mod- 



els in Figs. pHlQ (that do not occur in the corresponding 
Figs. ||-[|) grow because of the greater amount of attractive 
energy. The transition curves of Fig.|| and [l(] move to the 
left, those of Fig.|| to the right. 

If one of the component densities vanishes (let us for- 
mally set f^M := 0), then Figs.^ and 
formally reidentify fij = Qm {a\ = - 



10 



still hold. We can 



-1) and then the fig- 
ures represent the case of three component models with 
dust, radiation and one exotic fluid (which is described 
by VI2 and a^). As a difference to Figs. ||| and ^, now the 
densities Q,m and Q are the independent parameters. 
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Fig. 9. Worldmodels with dust, radiation and two A- like, 
exotic fluids {fix, a% = 0.9), {fl2, a\ = 0.5). In contrast 
to Fig.^ and corresponding to Fig.|| the effect of positive 
and negative energy densities have interchanged. For the 
solutions on the curves the same as for Fig. | holds. 
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Fig. 10. The domains of different solution types for one 
dustlike {fix, ax = —0.9) and one A-like fluid (f^, a.2 = 
0.9) as exotic additions to dust and radiation. The oscil- 
lating and recollapsing big-bounce models are separated 
by hill- type solutions (bold part of the transition curve) . 



7. Supernova constraints 

Employing standard least-square techniques we derived 
best-fit parameters for a number of models compared to 
the supernova data set of Perlmutter et al. (1999) (for de- 
tails sec Thomas 2000). We found that with a dust compo- 
nent Qm > representing 'ordinary' and cold dark matter 
in the universe the observations force a repulsive, A-likc 
component [fl, a) in our framework of description. Our 
best fit yields 

fl M = 0.86 ± 0.58 , fl = 1.38 ± 0.73 , a = 2 ± 1.9. 



According to Padmanabhan (1993) we write 

= 2.5 10- 5 /h 2 , with H = h 100^^ 

7 ' ' Mpc 

for the CMBR component. Using h = 0.65 one obtains 
fl 7 = 5.9 10" 5 . With this value the influence of the CMBR 
component is negligible. 

The great uncertainty in the value of the parameter 
of state a expresses that nearly any A-like, repulsive fluid 
can be fit to the data. When flM decreases then fl grows, 
while a decreases, too. For very low values of flyi , e.g. 
about flu ~ 0.03, close to the baryon density predicted 
by big-bang nucleosynthesis, one obtains a w 0.5. 

Among the models with a A-like repulsive fluid that 
fit to the data there are solutions without a big-bang be- 
ginning. These solutions lead to a ~ 0.5. A necessary con- 
dition for nucleosynthesis comparable to that of big-bang 
nucleosynthesis is to reach T ~ 10 9 K. Therefore one may 
restrict such world models to obey x m i n < 10 -9 . To ful- 
fill this condition flyi and fl have to be appropriately fine 
tuned which can be easily done. It is beyond the scope of 
the present work to check detailed nucleosynthesis scenar- 
ios for such exotic models, which would also depend on 
baryonic density and time scales at the nucleosynthesis 
epoch. 

A world model can only be considered as a successful 
fit if its age to exceeds the ages of the oldest known ob- 
jects. In Fig. [ll] we plotted thin lines of constant age in 
an OM-^-diagram for two-component models. The state 
parameters a of the exotic component with density fl are 
a = 2 (A-fluid) in the upper panel of the figure and a = 1 
(domain-wall like fluid) in the lower panel. The 90% confi- 
dence regions of the fits given in the diagrams are roughly 
encompassed by iso-age lines (for h = 0.65) correspond- 
ing to 12.5 and 15 Gyrs (for a = 2) or 15 and 20 Gyrs 
(for a = 1), reveiling the trend to larger to for decreasing 
a. When approaching the fat vertical-diagonal line from 
below, the world ages would increase towards infinity at 
this line, which separates big-bang models from bouncing 
models. 

Stars in old globular clusters are considered to be 
among the oldest objects whose age can be determined in 
a straightforward manner provided that calibration uncer- 
tainties (most importantly those of the RR Lyrae absolute 
magnitudes) are narrowed down. A recent revision of the 
ages of globular clusters yielded 11.5±1.3 Gyrs (Chaboyer 
et al. 1998). Allowing a formation time of globular cluster 
stars of at least ~ 1 Gyr leads to the condition to > 12.5 
Gyrs, which is fulfilled by the models inside the 90% error 
ellipse in Fig. [fl]. 

When taking into account recent results from the 
CMBR first acoustic peak which require a nearly flat uni- 
verse (de Bernardis et al. 2000) we find a tendency a — > 2. 
Consequently this favours the repulsive component to be 
the A-fluid. Our best fit with a flat model {k = 0) yields 

flu = 0.30 ± 0.04 , fl = 1 - fi M , a = 2± 0.04. 
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Fig. 11. This figure shows 68.3% and 90% confidence re- 
gions for world models that have been fit to the SN la data 
set from Perlmutter et al. (1999). The upper panel corre- 
sponds to 2-componcnt models consisting of dust (Om) 
and a A-fluid (density f2 and state parameter a = 2), 
whereas a composition of dust plus a domain- wall like fluid 
(density f2 and state parameter a = 1) has been utilized 
in the lower panel. The fat lines separate models of dif- 
ferent solution type, the thin lines represent models with 
constant age (in Gyr). 



f^M = 0.3 is in agreement with constraints from pecu- 
liar motion (Zehavi & Dekel 1999). Additional exotic fluids 
do not improve the goodness of fit and may therefore be 
considered as unnecessary. However, due to the increas- 
ing number of free parameters with each additional fluid 
one can fit models to the data that have for instance no 
big bang. Possible fits include models with x m i n < 10~ 9 , 
f^M = 0.3, f2 7 = 5.9 10 -5 that are flat and fit to the data 
with two exotic repulsive fluids, one of them A-like, the 
other dustlikc. 



8. Discussion and conclusions 

We have extended available classification schemes for CP 
world models by allowing for generalized ideal fluids with 
P = ()-l)f, 7 £ [0,2]. In particular, we have not re- 
stricted ourselves to only encompass incoherent matter 
('dust') and a A-field, though these are included as special 
cases. Other fields beyond A, which are repulsive due to 
sufficient negative pressure, are commonly applied in the- 
ories of the early universe (Kolb & Turner 1990). Since su- 
pernova la data require a currently acting repulsive com- 
ponent (e.g. Perlmutter ct al. 1999), whose physical nature 
is unknown, a more general consideration of possible so- 
lutions is in order. It has become popular to employ the 
quintessence ansatz, which implies a varying equation of 
state as compared to our approach (e.g. Caldwell et al. 
1998; Wang et al. 1999). This approach is less suitable for 
classification purposes. 

We here utilized a multifluid approach with fixed equa- 
tions of state (as defended by Turner & White 1997), but 
generously allowed for negative energy densities in addi- 
tion to the possibility of negative pressure. Such nega- 
tive energies need not necessarily be due to matter fields. 
They could, for instance, arise in the context of turbu- 
lence (cf. Marochnik et al. 1975a, b). From the beginning 
of relativistic cosmology negative energy densities were 
implicitely considered by allowing for the possibility of a 
negative cosmological constant A (see Sect. 2.1). 

General relativity for itself contains no strong con- 
straints on the contents of its energy-moment tensor. 
Hawking & Ellis (1974) summarized how singularity the- 
orems can be derived if one imposes certain 'reasonable' 
energy conditions. These energy conditions are not an- 
chored in deeper theory, but are rather assumptions that 
are apparently fulfilled by 'normal' macroscopic matter. 
Quantized matter fields can violate the energy conditions 
(e.g. Parker & Fulling 1973). The occasionally occuring 
'age problem' of CP world models, imposed by an uncom- 
fortably large Ho combined with a large age of the oldest 
known stars (often from globular cluster fits), leads eas- 
ily to a necessary violation of the energy conditions if one 
sticks to CP universes (Visser 1997). Part of the exotic 
components considered in our classification can also vio- 
late the weak and dominant energy conditions; the strong 
energy condition is violated by any repulsive component. 

As our analysis shows, negative energy densities as well 
as negative pressures do not lead to mathematical diffi- 
culties under the condition of homogenity and isotropy. 
The transition curves smoothly cross the borders of the 
quadrants. Instead they make four new types of solutions 
possible, that were not yet discussed. There are the 'hill- 
type', 'hollow- type', 'oscillating' and 'shifting-type' solu- 
tions (Fig. 1), all free of a big bang and therefore without 
a beginning at finite times in the past. 

In universes with up to two general-fluid components 
none of these new solution types appear except the 'oscil- 
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lating' type, and the relations in the density plane do not 
qualitatively differ from that of universes with dust and 
A. 

Essentially, fits to the luminosity distances as a func- 
tion of rcdshift to observed type-la supernovae suggest 
the presence of at least one repulsive cosmic component 
in addition to the 'ordinary' matter and radiation content 
(Perlmutter et al. 1999). We demonstrated that there are 
models with repulsive components which do not neces- 
sarily start with a big bang. However, if accepting the 
well-known successes of big-bang models, which include 
the explanation of the abundances of light elements and 
the CMBR, such bouncing models should be tuned to be 
successful in this regard as well. This requires at least 
that the scale factor falls short of x w 1CT 9 at the bounc- 
ing epoch to yield temperatures of ~ 10 9 K. This can be 
easily achieved with one additional exotic fluid. However, 
such a fit would lead to another 'fine-tuning' problem due 
to the lack of any known physical reason for the value of 
the tuning parameter. Since there is, moreover, no other 
reason for invoking such components we follow Occams ra- 
zor to confine ourselves on the minimal currently required 
number of cosmic fluids and do not pursue such possibili- 
ties further. It should nevertheless be emphasized that the 
so-called inevitability of big-bang models (e.g. Borner & 
Ehlcrs 1988) rests on such assumptions that the universe 
only consists of a minimal number of ft > fluids. 

With Qm and 7 as above and one repulsive compo- 
nent (O, a) we find that the universe will expand forever 
and has started with a big bang in almost all models that 
fit the supernova data. At the transition between models 
with big bang and those without, the age of the universe 
smoothly grows to infinity. An exotic O < component 
does not yield models in good agreement with the obser- 
vations. Our fits showed that instead a A-like fluid with 
O > seems very likely. In particular, if one relies addi- 
tionally on the first acoustic peak in the CMBR the uni- 
verse should be nearly flat. In this case the A-fluid itself 
fits best as repulsive component. None of the new solu- 
tions appears in this context. However, if assuming four 
components one may also tune bouncing flat solutions to 
supernova and nucleosynthesis conditions. 

Since the current Supernova la data reaching to z w 1 
are compatible with a large number of multifluid models 
further observations are needed to reduce the uncertainty 
in determining empirically the cosmic equation of state. 
The observational relation has to be improved by more 
accurate or many more observations (tightening the mean 
relation) and extended to higher rcdshifts where the model 
curves tend to diverge more. 
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